Notes on Quantum Mechanics 
Lectures by Prof. Barton Zwiebach 
MIT OCW Physics 8.05 


Herein find notes from Barton Zwiebach’s lectures on Quantum Mechanics, Physics 8.05 in MIT 
OpenCourseware. The full set of course notes is available on the MIT OCW web site, and 
they’re really complete and helpful. These here are just the main new ideas I’ve learned and my 
extrapolations from them. Zwiebach is an extraordinary teacher, and he’s clarified a whole 
bunch of concepts for me. 

Lecture 1: the Schrodinger equation 

At first, I thought Oh, no. Here’s a rerun of the wavefunction. My mind goes numb. 


This is different. Zwiebach is so enthusiastic and so clear it actually begins to make sense. 


Start with the Schrodinger general equation. 


ot 2m 0x2 


qe ere + V(x, 9) W(x,t) 


Note the symbols. Capital Psi refers specifically to the general equation. We’ll see later that 
little psi refers to the time-independent equation 


Ap(x) = Ep(x) 


Those parentheses on the right of the general equation by the way. That’s 4 , the time operator 
Hamiltonian. The total energy. 


Allowed potentials include, among the more common, square well (and related step functions), 
parabolic well, and delta function: 
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Develop the mathematical tools. I’ll drop the function parameters, for clarity. But keep in mind 
WY means V(x,t), while y means w(x). 


If Y represents a particle, we need maths to locate the particle in space and time. Y is 
complex. It has to be, given that i on the I.h.s. of the general equation. In order to locate 
particles in the real world we need real values. Logical choice is the metric from complex math, 
the density function 


p(x,t) = vv 
Consider position first. Main idea, as usual for unitarity in probability, is that the particle exists, 


with probability one, so we must be able to find it somewhere. In one dimension, the particle has 
to be somewhere along the real number line in the range minus infinity to plus infinity. So 


+00 
| Ye dx=1 


Note that in the integral, *W dx is the probability density of finding the particle somewhere in 
the interval (x,x + dx). 


All that is familiar. Unitary probability. Amplitude vs. probability. New is better understanding 
of the density function. I can sort of see it now; there it is on the Real line. 
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Figure. Probability density. Absolute value of area, orange, between a and b gives 
probability of finding the particle in that interval. 
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There are conditions at infinity. Both Psi and a have to go to zero at +00. Otherwise 


probability and momentum, among other things, blow up. 


Next is the continuity equation for the wavefunction. Here’s a puzzlement. It’s easy to visualize 
charge conservation, for example, in its continuity equation 


Any charge that escapes a region of space must have passed through the boundary of that region. 
In a one-dimensional system 
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Apply to the wavefunction. Consider an interval a to b on the real line. Any change in the 
density of the wavefunction in that interval must result from a density current. 


dp = J(a) —J(b) 


given the sign convention that J goes to the right. 


What is it that’s being conserved here? Conservation of probability. If amplitude translates in 
space — if the wavefunction is wiggling or a wave packet traveling — then the probability to find a 
particle must flow from one region to another. Probability of finding it somewhere remains 1, so 
the probability density has to increase in the neighboring interval by the amount it decreases in 
this interval right here. 


Next up, the operator H . In general, operators change the wavefunction. A is the time step 
operator. If the wavefunction is time dependent, e.g. if the wave is oscillating, then H updates 


the wavefunction to the next time step. 


Here’s the eye opener I mentioned at the beginning. Take the general Schrodinger equation. 


OW (x, t) (= a? 
ih ——— = | —-—~ + V(x, t) | @, ft) 


ot 2m 0x2 


Remove the time dependence, i.e. consider a wavefunction in a fixed potential and psi itself a 
function (remember we’ re talking mathematical functions) depending only on position as the 
independent variable. Rewrite. 


a. (=n a? 
ih W(x) = Im axe 1 V9) w(x) 


Parentheses on the right represent total energy. Relabel the differential operator on the left and, 
voila 


Hyp (x) = Ey(x) 
Note that Af is an operator and E is real. So this is an eigenvector / eigenvalue expression. We 
can solve the equation to find energy eigenvalues of the wavefunction, i.e. we can find the 


energy spectrum of a quantum system. Solve the differential equation. It is first order in time, 
so solution is pretty straightforward. 


ine) = EW) 


-iEt 
w(x,t)=e /n p(x) 
Neat! We can calculate how the wavefunction evolves over time. Draw the wavefunction. 
Cartoon animate it to watch it change over time. Update the cartoon frames by that exponential 


in energy. 


On to those energy eigenstates. The wavefunction will have particular energy solutions, b,W, 
which are basis states in a vector space. So a general wavefunction can be expressed as 


p= y bi; 
i=l 


where n is the number of eigenstates, i.e. energy solutions, for that particular system. The 
summation above is the spectrum of the wavefunction, e.g. the energy states of a hydrogen atom. 
Interesting physics occurs in degenerate states, when more than two eigenstates have the same 
energy. More on that later. 


The eigenstates are orthonormal, as expected in linear algebra. In mathspeak 


| Prin dx = brn 


All this is practically useful for calculating the general wavefunction 


is -iEj ey 
V(x, t) =) be hwy 
j=l 
and for calculating the coefficients of the eigenstates 


bea = | ax Wm W 


Pay attention what’s going on here. The summation above shows how to find the general 
wavefunction from the stationary function. All the time dependence is in the energy exponential. 
All the evolution, all the dynamics is in that energy function. The integral, finding coefficients, 
is all about orthogonality. Dot product (this really is just a dot product, in integral form) picks 
out the term in question. By orthogonality, all other products go to zero. 


Finally, expectation value. Given a general time-independent operator, A , what value can we 
expect on repeated / averaged measurements? Here ‘tis. 


(A)p(t) = i “wea dx 


Real value on the left. Functions on the right. So the integral is a functional, converting a 
function to a number. 
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Think about that a minute. Suppose you’re trying to find a particle’s position, looking for the 
expectation value of the position operator. Well, that argument in the integral is teasing out the 
likeliest position from the probability density, the product of those Psi’s. The integral is finding 
average value over an infinite range where total probability is one, so there’s no need for the 


usual 1 iy (b—a) coefficient out front. Particle is most likely to be found where the probability 


density is greatest. 


We’ll see later that we can interpret the integral as sum of the projections of the rotated state 
vector. The operator transforms the wavefunction. That’s what matrices = operators do. Rotate 
or stretch vectors. Assuming the wavefunction is normalized, then the dot product of the vector 
with its transformed self gives you the projection, how much of that wavefunction you can 
expect to find with that observation. Projections on vectors. Observables and how much you 
can expect to observe. 


Lecture 2: bound states 
Prof. Zwiebach starts out with theorems about bound states, i.e. states that go to zero at +oo . 
e They’re non-degenerate. No duplications of states at the same energy. 


e They’re real. 
e And they’re either even or odd functions. 
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The lecture includes corollaries and strategies for the proofs. See the course notes for details. 


Lecture 3: position and momentum 


Position and momentum are observables, i.e. they have a physical instantiation that we can 
measure. We can think of them also as different bases, different vector spaces describing a 
physical system. Prof. Zwiebach introduces the essential linear algebra. 


You can switch from one basis to the other using Fourier transforms. 


w(x) = | dpe '/n pp) 


—oo 


and 


Pp) = | dxe npn) 


—oo 


Note that the momentum operator acting on the x-basis wavefunction gives the associated 
eigenvalue relations, and vice versa for the p-basis wavefunction: 


a i = ipjx, 
[ ave ngo= ye nd) 
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as expected. It’s all in that exponential. 


One of the neat things I learned in this lecture is how to think of the wavefunction as a vector. 
Draw a one-dimensional w(x). Parse out the function over intervals € . The wavefunction has 
a value at each interval. Voila! A vector! 


y(0) 
I) 
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And the position operator is a matrix. Given 


kp (x) = x(x) 


Translate to linear algebra 
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Makes sense! And now I appreciate why the wavefunctions sit in such a huge (Hilbert) vector 
space! 


The rest of the lecture introduces the Stern-Gerlach experiment. Key is understanding magnetic 
moment and how a divergent external B field can separate spin-up from spin-down. 
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Figure. Stern-Gerlach apparatus. Collimated beam of ionized silver atoms traverses gradient of 
magnetic field, which separates spin up from spin down. Credit Prof. Barton Zwiebach, MIT 
OCW Physics 8.05. 
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Figure. Contrary to prediction of classical electromagnetism, electrons are only detected in one 
or the other of two states. Credit Prof. Barton Zwiebach, MIT OCW Physics 8.05. 
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Figure: Diagram of Stern-Gerlach results. Credit Prof. Barton Zwiebach, MIT OCW Physics 8.05 


11 


J 
ase | 
—r 
M> 
7 


a = 
ii _Sas v2 A Sx — 4 & SO% 
, 
~~ Zz : > , 7 
| = kh 


Figure: Series of Stern-Gerlach apparatus at different orientations. This is the heart of 
quantum mechanics! Credit Prof. Barton Zwiebach, MIT OCW Physics 8.05. 


See Zwiebach’s Lecture 3 for thorough discussion of spin one-half and the spin operators. I 
thought I had recorded those notes, but they’ ve disappeared somewhere. 


Lectures 4—7: Linear algebra 


I learned a whole bunch here. Zwiebach was away for a couple lectures, so Aram Harrow and 
William Detbold filled in. They zoomed through all the essentials. What a great review! 


They based their presentation on Axler’s text, Linear Algebra Done Right. I believe that title. 
Here are the main take-aways and what now makes sense that didn’t before (or that I just 
assumed I understood but really didn’t). 


It’s all about vector spaces and their properties. Components of vector spaces are fields and 
vectors, properties collected in CANNDII. Fields for our practical purposes are the reals and 
complex numbers. Vectors are of many sorts: polynomials, lists, etc. Addition of vectors 
commutes. Vector addition is associative. There is a null (zero) vector for multiplication. There 
is a negative (inverse) for addition. Multiplication is distributive. And identities exist for both 
addition and multiplication. 
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Hilbert space is a complex vector space that includes an inner product. Inner product is an 
operation that gives a field value when two vectors are multiplied; best example is the dot 
product. All physics occurs in Hilbert space. I think the physical implication here is that 
physical space requires a measure of distance, and that’s given by the inner product. 

A subspace U of vector space V contains the field and vectors u such that U C V and U 
itself is a vector space by all the definitions. Now it must be the case for some combination of 
U's that 


U,® U, ®...U, =V 
where the © mean the ‘direct sum.’ We’re adding subspaces to build the larger vector space. 


Now consider. Suppose each U comprises linearly independent basis vectors that span U . 
Then the direct sum spans V. The direct sum forms a basis for V. That’s the definition of a 
basis: a set of linearly independent vectors that spans the vector space. 


Operators are maps that transform vectors in a vector space V to other vectors in that same 
vector space. Operators can be represented by matrices, but matrices are more narrowly defined 
as operators in a given basis. The same matrix produces different results in different bases. 
Operators are functions and obey the distributive laws and linearity, but they do not necessarily 
commute. The commutation part is interesting; the commutator is a kind of eigenvalue relation. 


Suppose operator R = x and operator S = d/ gx: then, acting on the polynomial p = a 


[S,R]p = Ip 
I is like an eigenvalue here. Just like ih acts like an eigval in [x,p] = ih 
The dimension of a vector space 
dim V = dim null(T) + dim range(T) 


This is the Fundamental Theorem of vector spaces. null(T) are all those vectors that the 
operator T takes to zero (zero vector). range(T) are all the (surjective) transformations Tv 
otherwise filling the vector space V . Injective means an operator maps one-to-one in the vector 
space. Surjective means the map fills the whole vector space. 


Eigenvectors and eigenvalues of an operator are in the subspace of a vector space such that the 
operator acting on an eigenvector returns another vector in that subspace. 
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I’ve sure seen that before, but it makes more sense in the context of the spaces. Axler rules! 


Which leads to another useful theorem and perspective on the spaces. The eigen-subspace plus 
the subspace orthogonal to the eigens fills the next higher dimension. 


vV=U+Ut 


You can see that with the 2d x — y plane, spanned by e, and e, . Add the perpendicular e, 
and you’ve got R°. 


Maybe the biggest ‘aha’ of these lectures was Prof. Zwiebach’s explanation of Dirac’s notation. 
It derives from the inner product in complex space. The ket is a good ol’ regular vector. 
(Remember, whether a vector space is real or complex depends on the field, not the vectors.) 
The bra on the other hand (and this was the ‘aha’) is a map. It maps the ket to a complex value 
in the field. (Which is what happens with an inner product.) 


The bra’s compose an injective dual space to the kets. i.e. the bra is unique to the dual ket. 
Operations apply the rules for complex numbers. 


(a|b) = (bla) 
(a|P,b; + B2b2) = (Bb, + Bob2|a)* = (6, b,|a)* + (Bob2|a)* = By (al|b,) Pa Bz {a\b2) 


I skipped a step at the last, swapping bra and ket again. That’s why the complex star disappears 
on the last brackets. But you get the idea. 


Of course, you can represent the bra’s as row vectors; inner product then just becomes row times 
column vector multiplication. Handy! 


References are in the 8.05 notes, all you need to know, and of course in Axler. What great 
resources those are. 


Detbold and Harrow look like kids, but they are on the frontiers themselves. Harrow works in 
quantum computers and information theory. He’s a student of Chuang’s. Detbold runs 
numerical simulations of strong interactions, including conditions at the cores of neutron stars. 
Neat stuff! 
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Lecture 8: bra’s and ket’s and operators 
Now we get to the really powerful payoffs. Look at what the bra and ket logic can do. 


Take for example the definition of the dagger operator and its operations. A couple useful 
theorems for operators pop out. 


(Ttulv) = (ulTv) 
If we’re in a real vector space, then Tt = (T*)’. That’s one of the theorems. And you can use 
it both ways — find the 7’s if you’re in a real vector space or prove you’re working with a real 
vector space if you know the 7’s. Also note — just formalism — that (u|Tv) = (u|T|v) and, 
most important for practical purposes, we can treat the 7°s as matrices. 
The other handy theorem: 


(Ttulv) = (v|Tu)* 


Proofs for all these rest on index manipulation. I had it all clear in my head, but then it went 
fuzzy. 


Anyway, them’s important operators, the T’s. I think they’ll have a lot to do with projections 
and measurement, expectation values and such truck. 


Lecture 9: bra’s and kets (cont’d) 


Ha! Told you so. That there expression (u|Tv) = (u|T|v) is an expectation value! Zwiebach 
let it slip in passing. 


Main arguments in this lecture are extensions of the definition of a Hermitian operator. Tt is 
the adjoint, as defined by the bracket relations above, while Tt = T defines Hermitian adjoint. 
And for Hermitian operators all kinds of wonderful properties follow. For example, the 
eigenvalues of Hermitian operators must be real. So they can describe measurements. And the 
eigenvectors of different eigenvalues must be orthogonal. Proofs are pretty slick. 


Given T Hermitian and Tv = Av eigens. Start with bracket (v|Tv) and work both directions. 


A(v|v) = (v|Av) = (v|Tv) = (Ttv|v) = (Tv|v) = (Av|v) = A*(v|v) 
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Expressions on the two ends are equal, so the lambdas must be real for 2 = A*. Done! All that 
work finding Tt pays off! 


What about orthogonality of eigenvectors? Here ‘tis. Same kind of math manipulation. 
Given two different eigenvalues on the same Hermitian operator: 

Ay(V4|V2) = (A, |V2) = (TV |V2) = (V4|Tv2) = (¥4|AgV2) = A2(V1|V2) 
So (A; —A2)(v4|v2) = 0 


But we’ve stipulated the eigenvalues are different, so (v,|v2) = 0. The eigenvectors must be 
orthogonal! 


Onward to unitary operators. By definition a unitary operator is Hermitian and preserves the 
norm of any vector in the vector space. 


|Uv| = |v] 

From this, since U is Hermitian, 
|UtUy| = |Uv| = |v| 

So UtU =1. The identity. We’ll be using that. 
On to bra’s. Now here’s new. We’ve been used to linear algebra 

|av) = alv) 
Fine and dandy. Linearity on a vector space. But look what happens when we start talking about 
position and momentum, non-determinate variables. That is, position, a particle state, could be 
anywhere in an infinite dimensional Hilbert space. Proper interpretation of the bra, then, is that 
it represents the state of a particle, not the vector position of the particle. |x) represents a 
particle at position x. It is the state of the particle being at x. So 


lax) # al|x) 


Left side is a particle in the state of being at position ax. Right side a is the amplitude of a 
particle in the state of being at position x. There’s a difference! Similarly |X) is the state in 
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which a particle is sitting at 3-d vector position x. The bra itself is not the vector coordinate 
position. Subtle, eh? But make sense when you think about it. 

Lecture 10: Uncertainty 

Here’s the absolute coolest thing ever. Pythagoras rules! 


Question is: how do you figure the uncertainty in an operator? Suppose you’ re trying to find the 
value of an observable. Operate on the wavefunction. How certain are you of the result? 


Here’s how you find the uncertainty of a unitary operator A acting on a normalized 
wavefunction wW. 


Lpyyy> 
A Ayy> = (A- SAE NY 


A rotates |W) . First calculate the projection of Aly) onto the eigenvector |W). That will 
give the amplitude to find the transformed vector in that eigenstate. 


Projection. It’s an operator. It casts the shadow of a state vector onto a basis vector and thereby 
represents the component of that basis in that particular state. It’s a simple calculation. Using 
the |y) basis as an example, 


Pip) = Wl 


That’s it! Note that that the projection is a matrix, so a likely operator sure enough. And look 
here. Projection of A|w) onto the basis gives the expectation value in that basis. 
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IP) PlAlp) = |h)A) = (ADI) 


All that’s legal because the expectation value (A) is a real number (because the operator is 
Hermitian). 


Good enough. Expectation value gives the amplitude of the transformed vector on that basis. 


And, methinks, that there explains a lot! Look back at the Dirac definition of the wavefunction. 


p = | dx |x)(xhp) 


That’s the projection of psi onto the x-basis. Sum over all x and you’ve got the state vector = 
the wavefunction! Remember, x is a non-determinable basis. It’s the Hilbert space of all x. 
Each position itself serves as a basis. So the definition says the wavefunction is the sum of all 
the projections of the state vector onto all the bases. Same idea as a vector in three-space is the 
sum of its x,y,and z components. Makes sense! 


Projections onto basis states are components. Keep that in mind. 


But back to the uncertainty. Almost done. Look at the triangle. Because all the vectors are unit 
vectors, projection is cosine of the angle between A|w) and |y). Previous of Zwiebach’s 
calculations found that the uncertainty 


AAl) = Al) — (ADI hp) 


You have to have the J in there because we’re dealing with matrix operators. But take a look at 
the figure. AA|w) is the sine leg of that Pythagorean triangle. That’s the uncertainty, expressed 
as a vector. And Pythagoras immediately gives us the most useful relation for calculating 
uncertainties. 


(AA)* = (A’) — (A)? 


It’s all from the triangle. But the (A) isn’t obvious. Here’s the derivation. Remember the 
state vector is imaginary. So 


(Alp)? = (plA)(Alp)) = (lA? |p) = (A’) 
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Pretty neat! There’s a whole lot more algebra in the formal proof, but it’s all right there in that 
vector diagram. And a whole lot easier to remember. 


Note also the relation to statistics. That ‘mean of the square minus the square of the mean’ is the 
statistical variance. Variance. How widely dispersed are the elements in the sample. Applied to 
operators and states, I suppose it could represent the variation you’d expect if you performed a 
whole series of experiments on particles all prepared in the same initial state. 


Lecture 10: the uncertainty principle 


Variance and uncertainty, related by experiment. The uncertainty in repeatable experiments is 
the variance around the mean, the expectation value. Measure spin along the x-axis of zillions of 
electrons prepared with spin |z +) = |+) = |T) . We want to know the variance of our 
experimental results, our observations, our observable. What is the uncertainty in s, ? 


Think about that again. Prepare spin up along z. Measure along x with good ol’ Stern- 
Gerlach. Electron spin angular momentum is A] 2. Expectation value is zero (because it’s 


equally likely we’ll measure |x +) or |x —). We want to know how much variance there will be 
around the expectation value. Well, duh. Spin is quantized. There are only two possible 


outcomes to our measurements, + A] 2 or — A/ 2. Wariance is A] 2 - Uncertainty of the 


outcome is A] 2 - Uncertainty in the observable s, is Ay . Done. 


Let’s see if the math works. We want to find the variance in s, given state |z +). Use the 
formula. 


(Ass)? = (s,7) = (s;)? 


Convert to matrix representation to calculate those expectation values. 


x)= (Heit =f o15([° [=o 


(s,7) = (+|sy7|+) = [1 0] (*) F ; A = el 


Put it all together, 
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As expected! The logic works. Variance = uncertainty lower bound as expected from the 
quantization of spin. 


Okay. How about those famous uncertainty principles? AxAp = . and AFAt = ? Where 


did they come from? How can we understand them? And how can we figure out the mintmum 
of uncertainties, the minimum variance in our experimental results? 


Jumping right to the general formula, 


1 
AAAB > ((p| 14,8] lv) 


where A and B are Hermitian operators. Check that out with the known 


axAp > (| 5-[x.p1 bp) = («wl gin tw») = 3 


As expected. Note that the uncertainty is real. Now where did that general formula come from? 
Back to basics. We need to figure out that commutator in the expectation value. Begin with the 
Schwarz inequality. Product of norms of the two sides of a triangle is always greater than or 
equal to the product of the projection of a vector and its projector. 
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Schou Lneyvality 


U 
7} 
PU = |y<vlw> 


(av(v,v) 2 (uv) 
(uv)= VAY 


Apply that to our geometric interpretation of the uncertainty. 


Uncerdinty A . 
MW? 
nw? 


yyy) MP 
A Ayy>=(A- ADL 


AA) = Al) — (ADI hp) 


Add a second uncertainty to the mix. 
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AB|p) = Bl) — (B)I|p) 


Look at the geometry. 


Uncertainty Rancip |e 


IY> 
Ay (oB)* = Calad) 


The operators A and B are working in two different bases. They operate on the same state, 
ly) , but in different basis states. Position and momentum, say. Imagine the state B|y) is 
projecting out of the plane of the paper in that different basis, rotating its uncertainty AB|y) 
along with it. We’re interested in calculating (AA)(AB). We’ll use Schwarz. In preview, 
though, we’ll be interested the projection of AA onto AB . That tells all. 


We’re out to prove 


AAAB 


IV 


(| IA. B] bp) 
| 2i | 


Convert to the Schwarz form 


1 2 
aa*as? > ((pl 514.8] 1p) 
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Let 


Alp) — (A)I Ip) = fly) 


and 


Bl) — (BY hp) = gly) 


Schwarz gives 


(fIfglg) = Kf lg)? 


We focus on that rhs of the uncertainty inequality. First some basic, straight off the Argand 
plane complex geometry. 


((flg))? = Im(flg))*? + Re((flg))? 


Left side is complex, imaginary and real parts. Pythagoras separates them. Let’s see where that 
leads. 


i 
Im((fl9))* = 5 (fla) — (al f)) 


That’s just another manipulation of the complex. —(a + bi) — (a — bi) = b, the imaginary 


component of the complex z . 


Similarly, 


1 
Re((flg))° = 5 Uf lg) + (glf)) 


Now it turns out the real component usually doesn’t affect the Schwarz inequality. What matters 
is the imaginary component. That’s because the real term is always positive and always, well, 
real. It’s that differential in the imaginary component that’s tracking uncertainty. It’s the 
imaginary component that takes us out of the real plane and into the variance due to 
measurement in different bases. A measurement with meter sticks vs. a measurement with 
stopwatches, say. What’s the variance in outcome when you have to consider both the error in 
your meter stick and the error in the stopwatch. Measure a meter a bit too long and a time too 
short, then measure a meter a bit too short and a time too long and you’ve multiplied your 
variance. That’s what’s included out there in the angle theta between the uncertainty vectors. At 
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least I think so. Anyway, we’ll take Zwiebach at his word and focus on that imaginary 
component. Time to plug the operators back in. 


1 
Im((flg))° = 5 (pl (A — (4))(B — (BY) 1H) — Cpl (B — (BY) (A — (A) YY) 


Do the algebra. 


1 
5 HI — ADB — YI) — Cpl B— (BYE — (AY) 


1 
= 5) (Wl (AB — A(B) — (AB + (AXB))1)) — Cpl (BA — BYA) — (BDA + (B)(A)) hp) 


Good. Now take a look. Those are expectation values with mixed terms, scalars (other 
expectation values), and operators for which the bracketing psi’s will calculate expectation 


values. For example pull out that first A(B) term. We’re working in linear vector spaces after 


all. We can separate terms. 
(pl (—A(B)) Ip) = (BY pl (ADI) = —(B)*(A) 
It’s the product of two expectation values, itself a number. 


If you simplify all those manipulations inside the original brackets, 


1 
Im((f|g))* = 5 (plAB — (XB) 1b) — (op |BA — (B)A) hp)) 


By the rules, collect the terms in the brackets. 


1 1 
Im((f 1g)" = 5; (plAB — (AB) — BA + (B)A) |p) = 5 (pl [4, B] hb) 


As promised! Pretty slick! Look how that commutator is telling us the uncertainty. It’s the 
variance between measurement in different bases. 


Next up: derivation of the Energy-time uncertainty relation. Now we’ve got the general 
formula; all we have to do is plug in likely operators. 


Start with the time-dependent SE and a dummy operator, Q , that is time independent. 
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ih hp) = Hhp) 
Bracket the commutator of H with Q. 
(pI [H,Q1 bp) = Ql HQ — QH bp) = (1 HQ kp) — Cl QH Wp) 
= (HbL Q IW) (IQ Lp) = ~ihZ (LQ WW) — 1H WL Q kW) = —2in 1) 


Let’s check that. Work backward. 
7) _@ he 4) ol 
= (2) == WI ehh) = (Cal avy + Gb OL v) = (Hl @ Wb) — OHI IHD) 
1 1 
= FH wl HQ |W) — (pl QH |p)) = FH wl HQ — QH |)) 


= 5 (lH, Q1 IY) 
Looks promising, but I’ve lost some factors in there somewhere. 
Note that (Q) is time dependent. Plug back in to Schwarz. 
1 ° Ost 
(AH)?(AQ)? = (WI = 1H.@I Iw) =(-2=(0)) 


+ AHAQ = -n =a) 


Now identify a time interval 


_ AQ 
At= 30 
Ot 


~~ 


That’s the time interval required for the operator Q to change by amount AQ. 


Done! Well, factor of two or so. AHAt = h. Identify AH with its eigenvalue and we’ve got 
it. 
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AFAt=h 


Good enough. 


Lecture 11: Uncertainty and the spectral theorem 


Professor Z. checks out the position-momentum uncertainty with an example, 


2 


Pp 
H=— 4 
ia 


First convert to operators, then carry on. Gives pretty good prediction for the experimental 
values. (Not sure where he gets those experimental data.) See his lecture notes. 


On to diagonalization. If an operator has a complete set of eigenvectors on a vector space, then 


you can diagonalize the operator, eigenvalues along the diagonal with corresponding 
eigenvectors. Better yet, you can find an orthonormal diagonalization: eigenvalues along the 
diagonal and orthonormal eigenvectors. 

Here are the necessary transformations. Given 


Tv=Av 


with operator T € L(V) and v = {v;} aset of basis vectors on V . Then there exists an 
operator A such that Av; = u;, u = {v;} another set of basis vectors on V and 


Tu=dAu 


with T diagonal. That’s possible because you can choose v orthonormal so that 


Axx 
A 

Uj = AinV_ = 

Anx 

because 
1 0 
Vy = 0 V2 = 1 etc 

0 0 


26 


Better yet, you start with an orthonormal v and choose a unitary operator A such that u also 
is an orthonormal basis. That makes T as simple as possible: a diagonal matrix, all the 
eigenvalues along the diagonal, with orthonormal eigenvectors. 


How find that T ? Assume v orthonormal and A unitary as per above. 
Tv =aAv > TAv=AAv > A'TAv = A'AAv = dv 
A'Tu = Av 
There’s your diagonalized matrix, A’T . 


I’ve skipped over some subscripts. Those lambda’s, of course, are different in the different 
bases. But onward. 


Main idea in this lecture, and one of the key ideas for all Heisenberg-Dirac QM is the spectral 
theorem. If an operator T in vector space V has eigenvectors that form an orthonormal basis in 
V, then T must be a normal operator. By definition a normal operator has the property that 


(T',T] =0 


Normal operators include Hermitian and unitary operators. They’re the ones we need for 
quantum mechanics. 


The rest of the lecture involves proofs for the spectral theorem and various properties of 
diagonalization. I guess the practical takeaway for all this is if you can diagonalize an operator, 
then you immediately see it’s spectrum; all the eigenvalues sit there along the diagonal. And all 
the eigenvectors are just 1’s in the proper slot on the unit vectors. 


The other idea here is figuring out under what conditions you can simultaneously diagonalize 
two operators. Key requirement is that the two operators must have bases of the same dimension 
in the same vector space. Then you can relate one to the other by converting the bases. 
Diagonalize one and you know how to diagonalize the other by change of basis. 


Complications enter when you have a degenerate operator. Same energy, say, for several 
different eigenstates. You want to be able to find out what makes those states different if it’s not 
the energy. That’s the whole point of the simultaneous diagonalization of different operators. 
Represent one operator in terms of the other and you figure out what distinguishes them. 


pa | 


The matrix math is interesting. Turns out that the diagonalization of a degenerate operator gives 
you a block diagonal matrix. Each block represents the subspace of eigen’s with the same 
energy. That’s the conceptually important point. Degenerate eigenstates mean you have several 
vectors (states) composing a vector subspace all with the same energy. Each energy is a 
subspace of the larger vector space. Add up all those subspaces, a direct sum, and you get the 
entire vector space. Anyway, once you’ve figured out the blocks you can diagonalize the 
individual blocks with a new set of matrix transformations inherent to that subspace. There’s a 
whole lot of matrix calculations going on, but do-able. 


Lecture 12: Harmonic oscillator and quantum dynamics 

Harmonic oscillator. Again? For heaven’s sake, why? 

Well, we’ve got new tools and we can use them to pry further under the hood. Normal operators 
and the spectral theorem provide new understanding. Nuts and bolts. We get to see more of the 


blueprints, how quantum mechanics is put together. 


Start with the usual 


Factor. Purpose here is to produce a form H ~VtV so we can apply the spectral theorem. 


1 4) 4) h 
H =5mo*(2-—) (24+) 420 
2 mw mw 2 


Last term there, remember, results from the cross terms between factors and their commutation 
relations. 


Okay. So what? Well, those parentheses are normal operators! Let 


Normal operators. Their product gives a Hermitian operator, H. Their (orthonormal) 
eigenvectors provide a basis for their complex (Hilbert) vector space, so they are diagonalizable 
and their eigenvalues sit on the diagonals. We just have to read off the diagonals to find the 
spectrum of the harmonic oscillator. Can’t get any spiffier than that! 
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On to the nit and grit and an interesting aha! Define 


a: Pe pant ee PO ye 
Fal pe re ea) Oe 
Then 
1 
H = hw (ata +5) 


Tidy! Even tidier a number operator N = @'@. It’s going to tell us which energy state we’re 
in. It’s an integer. With that 


1 
H =hw (w a 5) 
2 
Some consequences of importance: 


[at,@] =—-1 and [@,a@‘] =1 


Check it out. With that, now here’s cool and illustrative. Remember working on the Energy- 
time uncertainty we found that 


ai 


Where Q is any general operator. Flip the other way, the commutator of the Hamiltonian on an 
operator is the time derivative of that operator with that quantum factor coefficient - . Well, 


looky here. 


[H,a] = —hwa 


You can derive that result from the essential properties of the ladder operators @, @t and their 
commutation relations. 
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Anyway, where have we seen that? 


—iwa 


S| & 
2 
I 


ewe 


2> 
II 


Plucked right out of 


w(t) = e-ty0) =e /n (0) 


It’s the phase operator that tells us how the wavefunction (state) evolves over time! Trumpets! 
Fanfare! @ is the operator that gives a state an energy kick up to the next energy state on the 
ladder. There’s the connection between state space and wave mechanics. 


Whew! I’m bogged down in the harmonic oscillator. Bogged down with the abstractions. 
What are the number operators? What’s |E)? What’s it mean |n) = = (a")"|0) ? What’s 


the relation of w(x) to all the brackets? Maybe back to the big picture to figure out all the 
parts. 


Take a look at that parabolic quantum well. 


Qhavitum harmon IC rcillute if 
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There’s the potential, built into the geometry of the well. There are the number operators, 
counting rungs up the ladder of the wavefunctions. There’s the lowest energy level, down at the 
bottom. There’s the momentum built into the wavenumber at each energy level. 

Focus on them counting operators. What information do they contain? Knowing the energy 


level (number) should provide enough information to fill in all the details. x , the only variable, 
is determined by the parabolic geometry. So the wavefunction, fitting waves into the parabola, is 


: : : a : 
determined, and momentum is also determined (by the ao operator). All you need to know is 


the energy level (number operator) and the potential. 


Good! See if we can make it work! 


We’ve got a number state |E). What’s the corresponding wavefunction? 
Figure out the expectation value (E|E). No help there, just a number. But here’s what we can 


do. We know that (atE|ate) = (N+ 1)hw (E|E) . If we just had a bottom line, the lowest 
energy level, we could bootstrap up to all the others. Well, we’ve got that. 


1 
Eo = (1, +5) hw 


Pe igae eo ik 
No = 0 by definition at the lowest energy state. So the ground state energy is 5 hw as per we- 


already-know-that. 


Now find the wavefunction at that lowest state from 


abo) = Jel. (2+) 


The lowering operator acting on the lowest energy state gives you no state at all. 


v.) =o 


Solve the differential equation and you’re set. You’ve got the ground level energy and you can 
construct the ladder up from there. 
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Here’s another perspective from R. Shankar’s text. Shankar dives into the weeds, but he’s 
careful to emphasize the essentials. Here are the essentials. 


@ and 4G" are defined from the Hamiltonian. From that, note [@, at] = 1, [A ; at] = at and 
[A ,a | = —@. That follows from the algebra of those operators. Note that the coefficients of the 
raising and lowering operators convert to units of energy. 


Most important of the operators for building the ladder: 
H=(ata+5)n =( +5)h 
=(@at+>)hw=|n+ 5 )hw 


where n is Shankar’s number operator. It records which rung on the energy ladder we’re at. 
Makes sense, @'a =n. Operate on an energy eigenstate with a takes you down a rung; then 
operate with @t and you’re back where you started, on the same energy level. Now climb the 
ladder. Energy of the ground state |0) is 


ieee o+;)a uly 
(n+5) hw =(0+5) ha => hw 


Energy of the second energy level 


and so on. 


: : : aH 
Next abstraction, rescale the energy into quantum units of hw. Operator H = or and 


eigenvalues € = o . Those are integer numbers. They scale the energy ladder of the harmonic 


oscillator into unit steps of hw. We can label the rungs of the energy ladder by those numbers. 
Call them |n) , energy states in the harmonic oscillator labeled by which rung of the ladder. 
That’s all we need to know to get the rest of the information about the state, e.g. it’s particular 
energy. That’s all we need to know to figure out the spectrum. 


At which we should take a moment to review Shankar’s take-aways from this study of the SHO 
in a parabolic potential. 


1. Energy is quantized. 
2. There are unit steps fw between energy levels. 
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ae 
3. Ground state energy is > , not zero. 


4. Wavefunctions leak out beyond the potential well, as per Adams’ discussion of 
exponential decay at the boundaries. 

5. Wavefunctions are either symmetric or antisymmetric in the potential well, i.e. cosine- 
like wavepackets or sine-like. 

6. Because of 5. above, the probability distribution for locating a particle does not follow 
the classical predictions. (Compare graphs of |nj) and |n,). Particle most probably 
at the center of the well in the ground state but probability zero to be found at the center 
in |n,). 


Ghavtum harmonic EciNutov 


Onward to the matrix representation of energy eigenstates in Hilbert space. 
Consider the Hamilton operator acting on a number state. 
Hale) = (€ — 1)ale) 
Work it out. 
Hale) = (aH + [A, al)le) = (ae — a)le) = (€ — 1)ale) 
as advertised. We’ve substituted the energy eigenvalue ¢ for the operator in the last steps. 
Take a close look. What that’s telling us is that (¢ — 1) is the eigenvalue for the operator H 


acting on the original state ale). Since the number operators identify the states, then ale) 
itself must be the same as the state |e — 1) up to a phase coefficient. Since the number operator 
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tells us everything we need to know about the energy state, if an operator transforms the original 
state aleé) into state (e —1)ale) then ale) must equal |e — 1) up to a phase factor. 
Hale) = (e — 1)ale) = Ale — 1) = (e — 1)|e — 1) up toa phase factor. 
ale) = C,|e — 1) 

Similarly we can show that 

Aatle) = (€+ Dale) = (€ + DCyaile + 1) 
The lowering and raising operators take us down and up the energy ladder. Up is okay. We can 
go forever up. Down must have a floor. A minimum ground state energy. What is it? Here’s 
one advantage of Shankar’s approach. 
At the ground state 


al€y) = 0 


Note the zero here is not a state, not an eigenvalue. It’s a representation of nada. Such a 
condition, a state below the ground state, does not exist. 


Well, if you act on a non-existing state with the raising operator, you still get nada. 
ataleo) =0 


Now apply the operator relations. 


atalea) = (A 5) leo) = (eo —5) leo) = 0 


So 
as 1 
Ey = 5 
or, translated back into units of hw 
Pi hw 
oD 
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oh : 
Ground state energy is = . That’s how we figure it out. Back and forth between number 


operators and the physical Hamiltonian. 
Now let’s go after those coefficients to the ladder operators. Here things get a little confusing 
with labels. Shankar (and Zwiebach, also) seem to be using the number labels and the energy 
level labels interchangeably. They both represent which step we are on the (quantized) energy 
spectrum. Anyway, using the number labels now: 
Given 
a|n) = C,|n — 1) 

it follows 

(nlataln) = (n — 1]C,"Cy|n — 1) 
By previous definition, @ta =n, the number operator. Rewrite. 


(n|n|[n) = (n — 1]C,'Cyln — 1) 


So 


Cy = Vn 


Coefficient for the eigenstate generated by the lowering operator is V7 , the square root of the 
energy level of the original state. Similarly, coefficient for the eigenstate after the raising 


operator is Vn +1. 
at|n) =vn+1 |n+1) 


And now — drum roll! — we can put it all together in operator matrix and normalized eigenstate 
vector form. Take a look. 


=) 


oo 
oO 
eo ets 
ae 
ooeo Oo 
cooNo Oo 
5 eee 
II 
+ 
8 


So when you raise the n = € = 3 eigenstate, for example 


mn IN 
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Oo NR 
NIW © 
Niu 


II 
~ 
S 
“Ni 
oo ogeceo 


NIL oO 


Pretty slick! Take a close look at how those row and column indices are working to raise from 
energy level 3 to 4. 


But that’s not quite the whole story. We’ve got to construct the spectrum from the ground state 
up, rung by rung up the energy ladder. And we have to normalize. Start with orthonormal 
eigenstates and you have to maintain orthonormal states. Not too bad, really. Here’s the grand 
finale. 


ee 
Ga 


|n) = 
Start with the ground state. Operate n times with the raising operator, as per the example 
above. Divide each time by the coefficient to maintain unit norm. Ta da! You’ve got the 
spectrum of the quantum harmonic oscillator! 


Best introduction to next idea is email I sent to Prof. Zwiebach re: his derivation of the 
Schrodinger equation from the time operator, U: 


I am working through your MIT OCW Physics 8.05. In Lecture 12, Dynamics, you derive the 
Schrodinger equation from the unitary time operator. It struck me that the result is (or sure looks 
like) a continuity equation, where the Hamiltonian operator is a 'current' and the norm of the state 
vector is a conserved 'charge.' 

nZ1y(oy = Ze ust we) 

ina WOES ap ol? 
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A continuity relation seems to make sense, a la Noether, with regard to conservation of 
energy. Anyway, I was curious if this notion has any merit. On superficial review of the 
literature, I find SE as continuity of probability, but I don't see any reference to this way of 
thinking about continuity of the time-dependent SE. 


I'm an old geezer, retired high school teacher, trying to figure out quantum gravity. I've been 
brushing up on QM. I've sure enjoyed your video lectures, and I've learned a whole lot. Thanks 
very much for making these ideas accessible to the rest of us. 


Here’s the full monty. Start with the Bloch sphere. 


0) 


Pole states: 


i+) = =5((0) + i|1)) 


|i) 75 (0) — i|1)) 


11) 


Figure. Bloch sphere. Credit Andreas Ketterer. 2016. Modular variables in quantum 
information. Thesis. 


A unitary time operator rotates the (normalized) state vector around the Bloch sphere in Hilbert 
space time step by step. Let U§ represent the operator that takes the state |w(0)) to |w(t)). 


Up (0) = hp(e)) 


38 


Time eperertov on the Blesh 2hbe 
Zz 


Figure. Time operator on the Bloch sphere. U updates the wavefunction in time increments. 


The unitary time operator is unique. It evolves any state to its next time increment. And by 
definition it’s reversible 


US! 
and 
uu, =t 
SO 
ut =u?" 


I’m using the indices here for general purposes. They could be any t1 and t2. 


Also 
UBUE = UR 


With that, we’re all set to derive the time-dependent Schrodinger equation. Start with 


a Be ac 
— pCO) = = (U5 1bO))) 


The only time dependence on the rhs is in the unitary operator, so 


= yo) = TO) 


We want the same eigenstate |y(t)) on both sides. 


= ee a 8 ys Ip(t)) = Spo) 


Call that operator on the rhs A. 


and 


Note that the lambda operators are currents! Just like 


0 veo 


yp") 


is the probability current. 


Claim is that the lambda operators anti-commute. 


A+ At=0 
Easy to show: 
T 
g t OU0 et Uo 
5; (UsUS") al = ap Ub t+ UG = At At =0 


Multiply both sides by a factor ih . That converts the lambdas to commuting operators. Then 
we're set. 


0 00S ac — 
AS |W) = tha — Uo Ip) = th ALP) 


Now ih A is a unitary, time-step Hermitian operator. What’s ina name? Call it A. 


0 Bs 
ih— |w(t)) = Abbe) 


Schrodinger! And cool thing is, by all appearances it’s a continuity equation. A hence Hf isa 
current. It’s a current in time. It’s the flow of time. 
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A OUG _ + 
H = ih—— UG 
ee 


That derivative-times-operator is a current just like the probability density current 


a 
[= “e w(x)* . The conserved ‘charge’ is the norm of the state vector. 


0 ss 
(POLin= WO) — A = 0 
A charge. A Noether current. Symmetry. It’s all right there. Energy is conserved. 


Well, hardly surprising. We’ve seen that from the commutation relations and a bunch of other 
things before. Still, pretty cool that it appears in Schrodinger. Or that you can derive 
Schrodinger from Noether. 


Lecture 13: Dynamics (cont’d) and the Heisenberg operator 


What’s next is to understand the unitary time evolution operator in terms of the Hamilton 
operator. Any calculation in quantum mechanics generally starts with a Hamiltonian. That’s the 
physics. From that we want to figure out U , to help with the maths. 


Finding the relation is straightforward. Return to Schrodinger. 


a a a 7 
ih lW(t)) = ih (US WW(O))) = Ale) = A (Us WO) 


Since the only time-dependence is in the U{'s we can write 


i Os, Rs 
ih aE U5 =H U5 
and solve! For a time-dependent Hamiltonian, 


-int 
UE=e Mh 


Similarly for a ‘slightly’ time-dependent Hamiltonian. 
sta to) 


trices 
U5 =e 


The whole-shebang Hamiltonian-time-dependence follows: 


4] 


Us = exp (-; [x it) 


That comes from the Taylor series of the exponential with the time ordering operator. See Prof. 
Z’s notes for details. Anyway, makes sense just by the looks of it. Time operator increments 
step by step over time, driven by the Hamiltonian. 


Onward to Heisenberg operators. By our previous definition, 
Ay = Ut" Agus 
H = Yo AsUo 


where the H and the S subscripts refer to Heisenberg and Schrodinger. Schrodinger operators 
are all the usuals, %,,H,etc. The dynamics now is shifted to the Heisenberg operator. 


By the definition, 
(P(CASHO) = (wCO)|UG*AsUE|(O)) = C(O) Auto) 


That’s handy. The Heisenberg operator allows us to choose an initial state to study the 
dynamics, and we can stick with that state through our calculations. 


This definition gives a bunch of handy relations linking Heisenberg to the usual Schrodinger 
operators. 


C5 =AsBs 9 Cy = US AUEUE B UE = UE ASIB SUS = UE" AsBUE = Cy 
and the commutator relations are unchanged. 
[An Bul = [US' ASUS, US" BsUS| = (US"AsUS) (US"BsUS) — (US BsUS) (Us "ASUS ) 
= (us'AsBsUs) = (US"BsAsUs) = AyBy — ByAq = [Ax Bu] 
sure enough! Note that out of laziness I’ve dropped the operator hat symbols. 
From this we can prove anew that in Ay = [Ay, Hy] . 


Start with Schrodinger. 


T 
, 0 Oat , (26 + OAs +, U5 


Convert the time derivatives of the unitary operators to their Hamiltonians. 
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- aA ze 
—ut' fA, Ut + ih (ué" = us) + Ut'A.A, Ut 


" e aA 
—ut" fi,Ag Ut + UST ASA, UE + ih (us" us) 


aA 
= [Ay, Hy] + ih (us" _ us) 


Now if Ags , the Schrodinger operator, has no time dependence, the last term disappears, and 
we’re left with the commutator relation as proved. 


None of this is real surprising, but it’s interesting to see how those Heisenberg operators work. 


Take a look now how the Heisenberg picture simplifies our understanding of the physics. 
Following comes from Prof. Z as well as Shankar’s text. 


First note the geometric relation between the Schrodinger and Heisenberg pictures. 
Schrodinger’s (vector) states rotate around Hilbert space against fixed coordinates (the 
eigenstates). Heisenberg, on the other hand, says the states are fixed and the coordinates rotate. 
It’s the basis vectors, as represented in the operators, that are changing over time. The physics in 
both pictures is the same. Calculations come out the same. Just a different way of looking at the 
world, and Shankar says there’s a whole lot of other models we might build. Lots of room for 
creative thinking. 
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Sibednger v6. Hesenbers aye RTOS 


Figure. Schrodinger vs. Heisenberg operators. Schrodinger rotates the wavefunction, 
Heisenberg rotates the axes. 


The Heisenberg picture simplifies our thinking about the physics. H makes the dynamics look 
the same as in classical physics. For example, in the quantum harmonic oscillator 


m= [e(-aamt 2") =[2(-a] 


Whoa! What happened to the second term in the Hamiltonian? Well, it commutes with X so 
disappears from the equation. 


oe = 1 [8p] = - 5 (pb + (2.81) 
dt h2m ~~ h2m PIP *) PAP 


where that last step pulled out one of the two p operators for each term; you have to calculate 
the commutator twice; once foreach p. So 


i he ha esa, Fos, ae eae 
a — Fp (& PIP + [%, p]p) = — Fp, (2thD) aa 
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Just as in classical mechanics! Heisenberg operators return the classical equation of motion! 


ao 
Same for © . 
at 


Same rationale here for dropping the first term in the Hamiltonian. It commutes with p . 


Op ahs i 1 21p al as i 1 2([p L|K + [p |X) 
AE AQ me Wr] =p omer (lp x1k + Lp, X)% 
Ue ae aA 
ap AD (2ihk) = —mw*k 


Shankar’s other observation deserves repeat. Back to the definition. Because the Heisenberg 
operators do all the lifting for time evolution we can solve the dynamics just based on some 
initial state, which presumably we can determine. Which is the whole point. We know a state to 
start with and we want to see how it evolves. 


Lecture 14: Coherent states 


Coherent states are replicates. Well, sort of. They share the same energy but differ in other 
observables. 


Take the quantum harmonic oscillator for example. Lowest eigenstate, the ground state, has 
energy = hw . But we could shift the apparatus a bit and the ground state over there is the same 


= hw as the ground state here. Coherent. 


Prof. Zwiebach starts the lecture with a review of dynamics: position and momentum operators 
as functions of time. 


X(t) = X(0)cos(wt) + ss (0)sin(wt) 
Mw 


B(t) = p(O)cos(wt) — mwxX(0)sin(wt) 


New addition is the Heisenberg dynamics of the ladder operators. 


ay = e ta 
and 
a," — piwt at 
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Essential tools for understanding coherent states are the translation operators. 


xo / 
ly =e h 


What it does is increment the position by an interval x) . Maybe a more consistent symbolic 
representation would be 


-ipxo 
T° Se hy 


x 


That is, the translation operator takes the state from position x to position x + Xo. 


Travslution operator 


Figure. Heisenberg translation operator on quantum harmonic oscillator. 


: Pee ; 
Note the relation to the momentum operator. p = -— . Translation by delta-x along the 


A 
quantum oscillator boosts the momentum up the ladder. Makes me wonder if the whole universe 
is a quantum harmonic oscillator... That seems consistent with the differential representation 


of the momentum operator anyway. 


Note the parallel to the unitary time operator. U§ takes the state from time 0 , some time we 
choose to call zero on our stopwatch, to some later time t. U§ is moving the state through time. 
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T° is moving the state through space. Which raises the question, is there some relativistic 


2 2 , 
relation T,°” — Uj” = S22 A metric? 


Anyway, back to the standard stuff. Given the symbolic representation, the algebra of the 
translation operator is pretty clear from the physics. 


SO 


That is, you’ve translated to a new position then right back to where you started. 


T,,Ty. =T 


x0“ Yo o+¥o 
That is, translation steps are additive. 


Note the representation for T;,, acting on the position operator. 


Think about that one. Talking operators here. lhs is an operator. It’s rotating basis axes relative 
to a state vector. 


(h(x) Tro! Teo WOOD) = WODIF WWD) + xo! 


Take a look at the state vectors. (See Figure above, translation on the quantum harmonic 
oscillator.) 


Ty1x) = |x + X9) 
Switching between state vectors and the wavefunction 
ly) > px) 
T,., |W) > W(x — x9) 


Note the minus sign. That’s the usual rule for translating functions across the coordinates. 
Minus sign if you move the function to the right. 


Now that we’ve got the tools, on to coherent states. In the quantum harmonic oscillator 


IZ) =T,,10) =e /n|0) 
0 
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That’s it. By definition, a coherent state is the ground state of the harmonic oscillator translated 


is position but still in the ground state. Slide the pendulum a bit to the right. Move the 
snowboard pipe a tad further around the hill. Energy states are unchanged. Coherent. 


Here’s an aside. How do you run a unit analysis quickly so you can understand, for example, the 


coefficients in the expectation values? It crossed my mind a good start is in the equivalents for 


the Planck constant. 
h= Et = px 
Play around a bit, makes sense. 


x 


standard units for kinetic energy, momentum times velocity. 


Or 


Hamilton equations of motion. 


Take a look then at the position and momentum equations of motion we’ve derived in our 
dynamics. 


X(t) = X(0)cos(wt) + ae (0)sin(wt) 
Mw 


— ae : : ; : ee 
That coefficient raed AU) should give units of x . Let’s see. First term on the rhs is fine. 
Second term needs some reckoning. Units. 


where I simplified at the second step using the harmonic oscillator V = 5 mw? x? and, third 


: 1 
step, the units w = = 


It works! No great surprise, but maybe it will help keep track of the coefficients. 
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Back to Zwiebach and the coherent states. A few more key ingredients. 
(XolXo) = 1 

General expectation values under translations: 

(oA) = (0Te,*AT.,|0) 
So, for example, expectation value for position under translation isn’t surprising: 

(Xo |X|Xo) = (O|(% + Xo)|0) = Xo 
But expectation value for momentum is a bit counter-intuitive: 
(XolB|Xo) = 0 


That’s because it’s a (Schrodinger) stationary state, wavefunction just sitting there. No 
momentum. Finally, expectation value for energy 


(oll) = (0|A|0) + 5ma®xq? = = hao +motx,? 


Energy in the coherent state is augmented by the potential at the (displaced) position x9 . 
Makes sense. 


And for future reference: 


~ 1a2le _ 22 
X09 |X“|Xo) = Xo + =—— 
(%o|2?|¥o) = Ko” + = — 


mah 
(Z0l67lZ0) = — 


(XolXP + PX|Xo) = 0 
Onward to the dynamics. We’ll use the good ol’ Heisenberg operators so we can access classical 
thinking. 


(p(t) |As|Xo (t)) - (%|Ay|%o) 
Try it out on the position operator. See what happens to the coherent state position over time. 


na 


X(0)cos(wt) + sin(wt) 


(Xo |Xy|Xo) = bi | = x,cos(wt) 
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As per the general results above, the momentum disappears. As it should. And looky! It’s the 
good old classical equation! The position oscillates around zero. 


Lecture 15: Coherent states and squeezed states 

Prof. Zwiebach calculates the general coherent state, which includes the ladder operators, and he 
explains the squeezed state. The math is complicated, and I won’t reproduce it here. Just lazy, I 
guess, but I’m getting antzy, want to head back to the frontier, quantum information and gravity. 


Time to get moving, finish up the QM review. 


Squeezed states are worth some discussion, though. I’ve always wondered what they were. 
Prof. Z. explains them well. Here’s the notion. 


Take a coherent state in the ground state of a Hamiltonian. It has an uncertainty 


where the subscripts identify the particular Hamiltonian. 


Now zap the system into a new Hamiltonian. Calculate the new uncertainty. 
h MW, h 
Ax2 = = |—— 
2M2W2 MzW 7 .|2M,W, 


M4W4 
= /— <1 
M2W2 


Now if 


as in if the energy of the second state is higher than the first, then the state has been squeezed. 
Think of a Gaussian. It isn’t as wide as it was to start with. It has a sharper peak. 
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Figure. Squeezed state illustrated as Gaussian function squeezed by boost to higher potential. 


Represent that transformation as an operator. If you want to build a squeezed state you need a 
squeezing operator. 


sy) = o~'/,(atat—aa) 


Note that it is quadratic in the ladder operators, and note the order of those operators in the 
exponent. Annihilation op’s have to be to the right, acting first on the state you’re squeezing. 
Otherwise the whole thing blows up, driving the state upward with successive creation operators. 
Now define the squeezed vacuum state 


|0,) = S(y)|0) 


Applications of squeezed states are really interesting. LIGO uses translation and squeezing 
operators to reduce noise in its detection system. The mirrors oscillate a bit because of thermal 
noise. That smears out the gravity wave signal. Solution: squeeze the detector photons so 
they’re less exposed to mirror fluctuations and translate them to where they should be if the 
mirror was absolutely quiet. 


la,v) = D(a@)S(y)|0) 


Pretty cool! Squeezed states do marvelous things. Perform sharper measurements. Send sharper 
signals. 
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Lecture 16: Photon coherent states and two-state systems 


Idea here is you can write the Hamiltonian for the electromagnetic field in a way that looks like 
the harmonic oscillator. 


1 
E= 5" + w*q’) 


Beware; E here is the electromagnetic field. Looks the same as the harmonic oscillator but 
without mass. It’s reasonable in units: [pq] = [h] . So convert to operators and declare 


where 
a= | ara 
q= ah a+a 
and 
, fhw.. ose 
p= =z (a a") 
0) 


1 1 
H = hw(ata+5) = hon +5) 

Z Z 
N now is the photon energy. Just like ladder operator stuff in the harmonic oscillator. What’s it 
all mean? Well consider the photon. It’s a quantum of the electromagnetic field. Think of it as 
a coherent state, a mode in the stupendous harmonic oscillator well of the electromagnetic 
universe. It has an associated momentum and potential energy, oscillating as it is between its 
potential boundaries. Mass on a spring but no mass. Just the spring oscillating. That said, we 
can assign the usual ladder operators to the Hamiltonian just as in the QHO. Same maths. 
All that said, we can think of the field itself as an operator. 


E = e,(e "a + eat) sin(kz) 
where the field is polarized along the Z-axis. 
Lectures 17 and 18: Two-state systems, ammonia and NVR 


Two-state systems include e.g. spin states and the ammonia molecule: you can capture them 
neatly witha 2 X 2 Hamiltonian matrix thusly: 


a2 


[90+ 93 Gi 


ie 
91792 Jor 


- | = Gol + 910% + 9295 + 9302 


In this mathematical structure, all the dynamics involves some kind of ‘precession.’ Magnetic 
moment in a magnetic field as the prime example, of course, but same maths describe the 
ammonia molecule and other two-state systems. 


Note the rubric to build models. 
1. Find a likely Hamiltonian 
2. Find the energy eigenstates and eigenvalues 
3. Find the expectation values 
4. Find the dynamics, i.e. the time evolution coefficients 


Ammonia is really interesting. It’s a two-state system, nitrogen either above or below the plane 


of hydrogen atoms, so it has an electric dipole. Put it in an electric field and you separate 
molecules by energy A above or below the ground state. 


Aronia molecale 


Figure: ammonia molecule. Nitrogen (green) oscillates across the plane of the three hydrogens. 
It is a dipole molecule with characteristic flip frequency. An electric field separates the two 
states. 


Eigenstates you can label, as per usual 
ick _ [0 
Iny=[5] and = [5 | 


Then the Hamiltonian becomes 


= Be 2 
L&E —A 


a2 


From there you can calculate dynamics, which (no surprise) includes terms like exp (=) and 


cos(wt) , where w is the Larmour (precession) frequency. From those dynamics you can 
calculate how long it takes, time 7, for an up state to flip down. And that, my friend, lets you 
build masers! 


Separate states with a gradient electric field. Send up state into a resonant cavity of just the right 
length such that the transit time = 7. That’s just right for the molecule to emit a photon of 
energy = 2A. Photons pile up in the cavity. Let them leak out and you’ve got a maser. Nobel 
prize for Townes et al in 1964. 


Pramonin MAIEr 


Figure: ammonia maser. An electric field splits ammonia beam into high- and low-energy 
states (relative to the field). High energy state enters resonant cavity with dimensions such 
that the beam drops to low energy and releases a photon as it traverses the cavity. 


It’s all right there in those matrix operators and state vectors and a little bit of math. (Well, quite 
a bit of math.) 


NMR uses the same maths tools. Apparatus has a twist to it, though, a rotating magnetic field. 
Put your target nucleus in a really strong, constant B, field. Adda rotating B field in the x-y 
plane. Nuclear spin precesses around B, and also around the (rotating) B, . Effect is to torque 
the spin axis down into the x — y plane. As it spirals down, it radiates at the frequency of the 
rotating B,. Tune the detectors to that frequency. You’re seeing mostly the hydrogens in water 
water molecules. Strength of the signal depends on the water concentration and the composition 
of neighboring molecules. You can get even more information from the damping time and 
relaxation time; how long does it take to spiral down, and how long to revert to alignment along 
Bes 
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Figure NMR. In aconstant external magnetic field B and rotating field in the x-y plane B, 
nuclear spin will precess from the z pole down into the x-y plane. The process emits cyclotron 
radiation as it drops into the plane, and that information is used to construct an image. 


Lecture 19: Tensor product and teleportation 


I finally get it! Tensor product is not really a multiplication. It’s a record-keeping system for 
multiparticle states. 


Main idea is that you can’t describe a multiparticle system just by listing the individual 

properties of all the component particles. It’s not enough to know the position and momentum of 
each individual particle. Those particles are correlated. Their wavefunctions interact. You have 
to keep track of all those correlations, all that entanglement. 


If V is the Hilbert space of one particle and W the Hilbert space of a second particle, then the 
Hilbert space of a system with both particles is V @ W. For example, given two spin-half 
particles, their vector spaces 


V=W={I+),I-)} 
and 


V @W = {|+) @I+), 1+) @|-),1-) @I+),1-) @|-)3 
By convention we’ll typically drop the @ between state vectors. 


VOW ={ lr) loll) 4 
=) 


Note that the dimension of the tensor product state is the product of the dimensions of the two 
component states. All the usual rules of linear algebra apply: scalar coefficients distribute and 
so do vector states. 


au®(v@w)=a(u@v) ®a(u@w) 


With those rules you can do all kinds of marvelous things. Like for example build a quantum 
teleportation system. The physical system in the following example uses spin states. Spin 
operators are unitary, so conserve probability and information. You implement the operators 
with varying magnetic fields. Them’s what goes into the Hamiltonians we call quantum logic 
gates. Them’s what’s the physical instantiation of the operators. The gates are magnets. 


Here’s an illustration of spin operators on Bell states. We’ll use them for teleportation. Define 
the Bell state 


1 
Ifo) = qa) + |-)I-)) 


Unitary. Normalized. Perfect. Now operate with the spin operators. Note that we have to use 
augmented operators, i.e. (/ ® 0,), since we have a two-particle system. 


a 
v2 


Think about that. The J in the operator preserves the state of the first particle. o, flips up to 
down and vice versa, acting on the second particle. Similarly 


lpi) = C@ %) @ Ibo) = =4)1-) + 1-)I+)) 


i 


Ip2) = U @ 2) @ Iho) Sn? se ese 3 2) 
1 
Ips) = U @ a3) @ Iho) = agi? Sy) 


Work backwards to the paired states. We’ll need those for teleportation. 


1 

I+)I+) = 5 (Ibo) + Ips)) 
1 

I+)I-) = 5 (Ib) — ile) 


1 
I-)+) = 5 (Ibs) + ilg2)) 
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1 
|-)I-) = 5 (lo) — I4)) 


OK. Teleportation. Alice and Bob share a Bell state 
1 
v2 


Alice grabs the state she wants to teleport to Bob. 


Ifo) = =(4)I+) + I-)I-)) 


Ip) = al+) + BI-) 


She interacts her states to form a tensor product |¢ ) ® |W) . The subscripts below track spins 
held by Alice and Bob and, C, the spins to be teleported. The tensor product represents the 
whole system: spins of all three particles, the entangled pair and the state to be teleported. 


1 
IPodaz ® - = yl alc 4)s + |+)al-)cl+)e ) 
+ RU Mal cle + |=jal=)el=)e) 


Well now. We can identify those leading A ®C states with Bell bases. 


1 
lDodas @ Wc = 5 Go) + |bs))aclt+)s + 61) — tlb2))aclt+)s ) 
1 
+5 B91) + ilb2))acl—)p + U0) — 163) acl—)e ) 


Regroup that last equation as factors of the Bell bases. 


1 1 
Ido)ap @ I)c = 5 |b) 4c (al+)e + B|—)p ) + 5 !b1)4c@l)e + B|+)p ) 


1 1 
+5 tbe) aclal—)e — Bl+)p) + 5 13)aclal+)e = B|=)a') 


Now look at that! Each term on the right is the AC Bell state times the associated spin operator 
(reflected in the signs; look close) on |W), . Alice has teleported |W) to B! All Bob has to do 
is operate on his state with the appropriate spin operator. Alice has to send him that information, 
which operator. Done! 


NR 


3 
IPo)az @ IW)c = > 16) ® Oj lW)p 
j=0 
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Be sure to check Nielsen’s quantum circuit for comparison. I think I finally see where that 
circuit comes from. 


On to EPR and the Bell inequality. Zwiebach presents the argument really nicely. I’ve been 
confused about what is local realism and other such truck. Turns out it’s right there in the two 
assumptions Einstein insists on: 


1. Any measurement reflects a reality of the system. i.e., if your measurement determines 
that a particle has spin up, then that particle most definitely had spin up before the 
measurement. QM, of course, says that the particle was in a superposition of states 
before the measurement. 

2. Conditions far away cannot affect measurements right here in the lab. QM, on the other 
hand, says particles can be entangled, i.e. correlated, over vast distances. 


Bell’s inequality established what’s what. It’s straightforward in Zwiebach’s presentation. Key 
is that for a spin system the probability of measuring that spin lies along some axis at angle theta 
from the reference axis 


Check the math. Maybe it’s in Adams’ notes a while back. But I’m pretty sure I verified this 
myself. 


OK. Here’s Professor Z’s argument. Consider an experimental apparatus that can measure 
particle spin along any of three axes. Prepare entangled pairs. The table below lists all the 
possible entangled states. + and —’'s are spins along the three axes a,b,c . Columns list the 
measurement outcomes. State labels are arbitrary, just a counting device. 


state particle A particle B 
N, +++ --- 
N2 ++- —-+ 
Nz +—+ -+- 
Ng +—- —++ 
Ns —++ +—- 
No —+- +—+ 
Nz —-+ ++- 
Ne —-- +++ 


Figure the classical probabilities, what EPR predicted assuming local realism. Calculations list 
spin state for particle A followed by state for particle B. 


a(+)b(+) = N3 + Ny 


a(+)c(+) =N,+ Ny, 
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b(+)c(+) = N2 + Ng 
From those equations, it’s clear that 
Pla(+)c(+)] < Pla(+)b(+)] + P[bG)cG)] 


That’s the classical prediction according to EPR. But QM says it ain’t so! Suppose the angles 
between axes are small and the b axis lies between a and c. 


Bell tects 


——_  ——o— 


Bob 
S}—--~ 10 
Abe G 


() 


> 
ra 
c 


Figure. Experimental test of the Bell inequality. Alice and Bob independently measure spin 
orientation of their particle from an entangled pair. Bell inequality is obviously violated at small 
angle differences between the two measuring apparatus. 


Those probabilities are 


Not so! At small angles 
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That’s the QM prediction. Alain Aspect and many others have carried out the experiments. QM 
wins. 


Lectures 21-23: Angular momentum 


There’s a whole lot of really dense definitions and proofs here. Main point is to develop angular 
momentum operators, general J's that include all the maths of orbital angular momentum, spin, 
and everything. 


Main point, after a lot of work (which I really need to figure out sometime) is 
iJ; |] = ih €ij Sn 


Quantized. And note that A has units of angular momentum. Like x-p and E-t. Think 


about those relations. Make them all operators. ¥ + is obvious; that’s the classical angular 
momentum. £-é needs some more thinking. 


Anyway, after all the maths gymnastics we end up with a Hamiltonian for a spin in a central 
potential. Like an electron in the electromagnetic potential of its nucleus. 


nh? 1 a” h? 1 O oe: 1 @? 
— sind — + V(r) 


~ 2m 7 or?’ 2mr2\sind 06 00 + Sin? 6 6 ah 


Note all those accelerations along the radial and angular directions. That seems a quick 
shorthand to think about it, anyway. 


From that Hamiltonian you can show that the energy is quantized with quantum numbers j in 
increments of sh components along the z-axis and total momentum between +/j. I think that 
about captures it. 


The generalized wavefunction, with all the angular and radial terms collected into Y and u 
operators (via algebra to combine all the messy coefficients) 


Wem = ta) : —— + Yim (8, p) 


p isa function of radial distance. The wavefunction depends on radial distance and spherical 
angle. There’s just a whole lot of calculation goes into figuring those. From that we get the 
Hamiltonian 


he 07 


Im ape aay Ue(T) + VeppUg, = Eup, 


where the effective potential includes a term for centrifugal force 
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h71(1 + 1) 
Vers VO ae 
Take-home from all this — Zwiebach’s words — is summarized in the graph for orbital angular 
momentum. Note we’re talking orbital L here, not spin, so energy levels are unit quantized and 
not half-integer. Beyond the lj) (i** energy level with zero angular momentum, as in 

st,s*,s3, etc. orbitals) the states are degenerate. So, for example, there are three states at each 
l;, , one for each of the three 1, and with components labeled by m in the algebra. (Note that 
there are m = 21+ 1 z-components of angular momentum. See why?) 


Brgy spectros for oferta Grgluy momettun 
bys 


anelaner 


\ 
lev ‘ 


Take a look here how the spectrum is higgledy-piggledy, not nice and neat like the QHO 
spectrum below or the Hydrogen spectrum. 


Interesting also is that calculating the wavefunctions on a 2-d surface requires 3-d angular 
momentum. In fact, the L,,, operators emerge naturally from the algebra. (Don’t ask me to 
demonstrate that right off. Check out the lecture notes.) 


What do those wavefunctions look like? 3-d angular momentum? Well, they’re Bessel 


functions. Those are the 3-d standing waves, e.g. representing the modes of vibration of our Mr. 
Sun. Below is a vibrating membrane, but you get the idea. 
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Figure: Vibrating drum Bessel function. Credit Slawomir Bialy. Bessel Function. Wikipedia. 


Turns out wavefunctions in a uniform spherical potential well are a mess, no pattern. But the 
spectrum of a 3d quantum harmonic oscillator is nice and tidy. Here’s the algebra. Hamiltonian 
has the same form expressed in 3d. Then it’s all numerology. 


a 
2m 


H= 


We build the spectrum of the 3d QHO like we did the Id oscillator but with two more sets of 
ladder operators. And, note this, those operators form entangled states. The states in the QHO 
are tensor products. That’s interesting; before we were entangling particles. Now we’re 
entangling operators. 


Start at the ground state, no angular momentum. |wW) = |0). In 3d that has energy = hw . Now 


spin it up to one unit of angular momentum. There are three possible spin states, al |0), 
‘ : 5 
al |0), and a!|0) all with 1,. Degenerate in energy 5 hw . 


l, gets trickier. Six possible states, combinations of the creation operators adding to two units of 
half-spin. 


atat|0), atal |0), atah|o), ata}|o), atat|0), ata}|0) 
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Well now. There’s a problem. Six states totals energy “hw . Ain’t no such energy on the 
spectrum. We’ve got to split that degeneracy. Solution is five 1, states and one [po all at 
E= “hw . Onward and you get a spectrum that looks like 


Specttam of 3d std 
cle te thas 


S 
1S; 2a 
ip 2 
ly, T ys 
Y, )O0 
4, + 6 

3 
7), — | | 


But hold on here. We built those six states in 1, . How did we end up with a state in ly ? The 
answer (I think — Prof. Z. didn’t address this directly) is entanglement. Entangle all those l, 
states and you get an isotropic system among them. It’s spherically symmetric, no preferred 
direction. 


On then to the hydrogen spectrum. With all these tools available it’s simple! I’d assumed H 
would be a colossal maths challenge. Not so! It all spills out of the Hamiltonian of a central 
electric potential. 


a2 2 h2 2 
2m f 2mr2 F 


We can solve for the Bohr radius immediately. Set the potential and kinetic energies equal. 
Solve for radial distance. 
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h2 
ag = — — 
2me2 


Use that to solve the ground state energy: plug ap into the potential, then ladder up the energy 
spectrum with the angular momentum operators! Same drill. 


e? 1 


a 
mu 2a, n2 


Surprising is how neat and tidy is the hydrogen spectrum. All kinds of degeneracy. You can see 
all the orbitals right there in the spectrum. 


Hydegen Spectum 
cle) He S 


7, 2y 
IY >| 
\, = — TIS 
bp — i ean 
= +6 
HY 
+ 3 
i— 


Lecture 24: Intro to perturbation theory 

Here’s interesting: the Feynman-Hellman theory. RPF figured it out as an undergrad. 

Idea is that if you know the state of a system and you tickle it, say, with an extra potential you 
can add the perturbation to the operators on the initial state and get a close approximation to the 
perturbed state. Zwiebach’s example here is fine splitting in the hydrogen spectrum due to the 


magnetic moment of the electron. 


Feynman-Hellman says 
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OE(A) dH(A) 
= A A 
= -{v aloe v@) 
where A is the perturbation of the magnetic moment. Simple enough, and it makes sense. You 
can check it out with the bra and ket algebra. 

Now suppose Hey = Ho + a“ dA =H,+AHy . Inthe case of fine splitting, all we have to 


calculate is that second term. It’s the change in energy due to interaction of spin orbital 
momentum with the magnetic moment of the electron. That’s what we’ll add as a perturbation. 


2 


h 


where the orbital L and spin operators can be calculated from our previous calculations of the 
orbital radius and the fine structure constant. Not bad! 

Lecture 24: Spin-orbit coupling 

Question is: what happens when the system includes multiple components. For example, the 
hydrogen atom has a central potential, the electromagnetic field anchored on the proton, plus the 
magnetic moment of the electron in that field. The motion of the electron induces a magnetic 
field, and fu, interacts with that field. That’s the perturbation conditions introduced above. 


Now take a look at the states and their energies. 


Represent the possible |/,m) orbital angular momentum states. / is total angular momentum; m 
is the Z component. There are three possible states for 1 = 1. 


1,1), |1, 0), and |1, -1) 


And there are two possible electron spin states. 
1 ) d FE 5] 
re) aaa oe 
Altogether, then, there are six possible spin-orbital terms in the Hamiltonian. 
111) [5.5] 
: 202 
1,1) @ [,-3) 
; 2’ 2 
11,0) [5,5] 
; 22 
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1,0 

| f ,® E»-5) 
11 

|1,-1) ® 5.5) 

1,-1) 

| > ® -. — 


The top and bottom states have energy = hw . The middle four form a degenerate multiplet with 


energy = hw : 


It’s all wavefunctions, but we can see the picture in a cartoon. 


Spin -bit angalar Momecton 


i 


Figure. Spin-orbit angular momentum. Depending on magnetic moments of the nucleus and 
electron, relative orientation of spin and orbital angular momentum split the spectrum into fine 
and hyperfine spectra. 


Lecture 26: The hydrogen spectrum 


Done! Last of the lectures! And what a great lecture series! Prof. Zwiebach and his assistants 
do a marvelous job presenting the quantum mechanics. 
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Can I do the calculations? No, not well. I need to really sit down and practice. Do the problem 
sets. Do the exams. 


But I think I have a whole lot better understanding of the concepts. I learned tons about vector 
spaces and operators. I learned tons about wavefunctions and about Hamiltonians and complete 
sets of commuting operators. Tons. 


The hydrogen spectrum requires a bunch of vector algebra. Blackboard after blackboard of 
equations. What it all comes down to, I think is: 


1. Look at those 3d Bessel functions. In them you can see the angular momentum, total and 
z-component. 
2. Include the additional angular momentum from the Runge-Lenz vector. 


Runge-Lenz R points along the major axis of an ellipse. It has constant magnitude, depending 
on p. If the orbital is precessing, as in the presence of the orbital B field, then that precession, 
captured by Runge-Lenz, contributes to total angular momentum. 


Suppose the electron has an orbital angular momentum L and spin S$ with total f=L+S . 
Now suppose the whole system is precessing, major axis of the atom revolving around the 
nucleus. We have to include those effects, the precession, in calculations of angular momentum. 
That takes the form R x f. There’s a whole bunch of commutation relations in there, and it all 
shows up in the spectrum. 


Meantime I’m off to do taxes and get ready for a Grand Canyon trip. Hasta luego! 
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Epilogue 


Grand Canyon has been postponed because of CoVid-19. School and community life, all have 
been postponed. 


But there’s still quantum mechanics to learn. I’m embarking on Scott Aaronson’s lectures for 
quantum computing. Hasta hasta luego! 
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